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Dependence of CMI Growth Rates on Electron Velocity 
Distributions and Perturbation by Solitary Waves 

R. L. Mutel,^ W. M. Peterson,^ T. R. Jaeger,^ J. D. Scudder,^ 

Abstract. We calculate growth rates and corresponding gains for RX and LO mode 
radiation associated with the cyclotron maser instability for parameterized horseshoe elec- 
tron velocity distributions. The velocity distribution function was modeled to closely fit 
the electron distribution functions observed in the auroral cavity. We systematically var- 
ied the model parameters as well as the propagation direction to study the dependence 
of growth rates on model parameters. The growth rate depends strongly on loss cone 
opening angle, which must be less than 90° for significant CMI growth. The growth rate 
is sharply peaked for perpendicular radiation (fcy = 0), with a full- width at half-maximum 
1.7°, in good agreement with observed k- vector orientations and numerical simulations. 
The fractional bandwidth varied between 10""' and 10~^, depending most strongly on 
propagation direction. This range encompasses nearly all observed fractional AKR burst 
bandwidths. We find excellent agreement between the computed RX mode emergent in- 
tensities and observed AKR intensities assuming convective growth length Lc ~20-40 
km and group speed 0.15c. The only computed LO mode growth rates compatible ob- 
served LO mode radiation levels occurred for number densities more than 100 times the 
average energetic electron densities measured in auroral cavities. This implies that LO 
mode radiation is not produced directly by the CMI mechanism but more likely results 
from mode conversion of RX mode radiation. We find that perturbation of the model 
velocity distribution by large ion solitary waves (ion holes) can enhance the growth rate 
by a factor of 2-4. This will result in a gain enhancement more than 40 dB depending 
on the convective growth length within the structure. Similar enhancements may be caused 
by strong EMIC waves. 



1. Introduction 

The electron cyclotron maser instability (CMI) is widely 
recognized as the generator of auroral kilometric radia- 
tion (AKR) from the Earth's magnetosphere (see Treumann 
(2006) for a recent review). It is also been invoked to explain 

radiation from many other astrophysical sources, includ- 
ing Jovian magnetospheres [Zarka, 1998], microwave bursts 
from certain solar bursts [Melrose and Dulk, 1982; Dulk and 
Winglee, 1987], flare stars [Bastian and Bookbinder, 1987], 
white-dwarf binaries [Willes and Wu, 2004], and even cx- 
tragalactic radio jets [Begelman et al, 2005]. The free en- 
ergy of the maser derives from mildly energetic downward- 
directed electron beams which are accelerated by parallel 
electric fields and pitch-angle scattered by magnetic field 
gradients at lower altitude. The resulting distribution func- 
tion, sometimes referred to as a 'shell' or 'horseshoe' distri- 
bution, has a large crescent-shaped shoulder of positive ve- 
locity gradient {df/dv± > 0), resulting in strong non-linear 
growth of waves near the electron cyclotron frequency. The 
dominant free-space emission is the right-hand circularly po- 
larized extraordinary (R-X) mode, whose spectral power has 
been measured as large as 10"'* m~^ Hz"'^ [Ergun et al., 
1998a] for bursts with typical bandwidth of one KHz. This 
is equivalent to an energy fiux 2- 10~* W m~^ , which is ~1% 
of the kinetic power of the electron beam, a very high effi- 
ciency since many of the beam electrons precipitate to the 
ionosphere and do not contribute to the CMI instability. 
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Early papers which identified AKR emission with a CMI 
mechanism often assumed a loss-cone electron velocity dis- 
tribution (e.g. Wu and Lee [1979]; Lee and Wu [1980]; Wong 
et al. [1982]; Omidi and Gumett [1982]; Omidi, Wu, and 
Gumett [1984]). A key insight was that the relatively small 
relativistic mass correction applied to the cyclotron reso- 
nance condition resulted in a resonance ellipse rather that a 
straight line in velocity space [Wu and Lee, 1979]. The ad- 
vantage of the ellipse is that it could intercept a significant 
fraction of the 'inverted' {df /dv± > 0) regions of the loss- 
cone. However, Ormdi and Gumett [1984], using S3-3 elec- 
tron distributions and ray-tracing models of efi'ective path 
length, showed that the computed path-integrated growths 
were too low to account for observed AKR intensities. An- 
other problem with loss-cone models is that in order for the 
resonant ellipse to intercept the 'shoulders' of the loss cone 
(where df/dv± > 0), the ellipse must be significantly offset 
from the origin, which in turn requires the k-vector to have 
a significant upward directed parallel component. This is in 
conflict with in situ AKR observations in the auroral cav- 
ity, which indicate that the k-vector is perpendicular to the 
magnetic flcld direction within a few degrees [Ergun et at, 
1998a]. The loss-cone model became even more problematic 
when Louarn et al. [1990] reported that Viking observations 
of electron distributions in the auroral cavity having loss- 
cones with large perpendicular velocity gradients were often 
not associated with any AKR emission. 

A second key theoretical insight was that relativistic cor- 
rections also play a crucial role in the dispersion relation, 
even for mildly relativistic electrons [Pritchett, 1984a, b; 
Strangeway, 1986]. In the low-density environment of the au- 
roral cavity, the relativistic correction causes the RX-mode 
cutoff frequency to be shifted below the cyclotron frequency. 
For nearly normal propagation, this produces a new rela- 
tivistic mode which is highly unstable [Pritchett, 1984b] . It 
is this mode which dominates CMI growth. 
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In the mid-1980's, detailed simulations of the CMI mech- 
anism lent strong support to the idea that the growth rate 
was driven by an electron distribution function with a shell 
or horseshoe geometry rather than a loss-cone [Pritchett 
and Strangeway, 1985; Winglee and Pritchett , 1986]. This 
was confirmed in the early 1990's when in situ measure- 
ments indicated that the distribution function in auroral 
cavities does indeed often resemble a horseshoe shape [Roux 
et al., 1993; Delory et al., 1998; Ergun et al., 1998a; Strange- 
way et al, 2001]. The horseshoe shape is a consequence of 
Earthward-directed beams of 1-10 keV electrons which are 
accelerated by parallel electric fields and pitch-angle scat- 
tered by the converging magnetic field. The origin of the 
parallel electric field is controversial, but may be a result of 
mid-cavity double layers [Ergun et al., 2004]. The horseshoe 
distribution results in much more robust CMI growth since 
the resonant ellipse can be centered on the origin (fcy ~ 0) 
and hence can intersect a much larger volume of velocity 
space with a positive velocity gradient. This realization has 
resulted in a number of papers calculating CMI growth rates 
using specific horseshoe distributions [Bingham and Cairns, 
2000; Pritchett et al, 2002; Mutel et al., 2006]. 

Both left-hand ordinary circularly polarized (L-0 mode) 
and R-X mode emission arc detected by remote spacecraft. 
Internal Z mode is also detected in the source region [Gur- 
nett et al., 1983], but may result from conversion of RX- 
mode radiation propagating normal to the auroral cavity 
[Pritchett et al., 2002]. Although the R-X mode almost al- 
ways dominates observed AKR emission, there is often sig- 
nificant power in the LO mode, usually in the range 1% - 
10% of the RX power [Mellott et al, 1984; Benson et al, 
1988; Hanasz et al., 2003]. When L-O and R-X mode are 
seen together, their observed fluxes are often highly cor- 
related, with L-0 mode flux densities about 2% of RX- 
mode [Mellott et al, 1984]. This suggests that LO-mode 
could result from mode conversion of R-X mode radiation, 
perhaps by reficction at sharp density boundaries [Calvert, 
1982; Hayes and Melrose, 1986]. More recently, Hanasz et al 
[2003] reported that L-0 mode AKR is detected much more 
often at low magnetic latitudes (\m ^ 60°) in the afternoon 
auroral sector compared with the evening sector. They at- 
tributed this to a propagation effect: R-X mode originates 
in the evening sector where AKR bursts arc most often lo- 
cated [e.g., Huffet al, 1988; Mutel et al, 2003, 2004]. How- 
ever, some fraction of RX mode radiation is mode converted 
to LO mode (e.g. by refiection from density boundaries), 
and propagates along largely unrefracted ray paths over the 
polar cap, allowing some spacecraft locations to intercept 
only L-O mode ray paths. This idea was first suggested by 
Hashimoto [1984] and was also soon in DE-1 spacecraft ob- 
servations of R-X and L-0 mode AKR [Mellott et al, 1984]. 

The question of whether L-O mode AKR is primarily 
generated directly by the CMI mechanism, or whether it 
is mode-converted from RX-modc, remains an unresolved 
issue. Since the CMI growth rate is proportional to number 
density, in regions of relatively high density, L-0 mode is 
favored since the R-X mode is strongly suppressed by even 
modest Wpe/nce ratios. In the limiting case of no cold elec- 
trons, the relativistic RX-mode cutoff frequency is 



+ 



1/2 



(1) 



where 7 is the Lorentz factor of energetic electrons. This 
cutoff exceeds the maximum allowable frequency for CMI 
growth (fice) at very modest densities. For example, for 
5 keV electrons (7 = 1.01), CMI generation of RX-mode 
is suppressed for Wpe/f^ce > 0.1, corresponding to a lim- 
iting density = 3 cm~^ at frequencies below 225 kHz. 
Since the cutoff frequency for LO-mode is uipe <C fice in the 



lower magnctosphcro, it can contirme to grow linearly with 
increasing density if the CMI resonance ellipse (whose topol- 
ogy depends in part on density) intercepts a suitable region 
of the velocity distribution. 

AKR exhibits a bewildering variety of fine structure in 
dynamic (frequency-time) spectra. Curnett and Anderson 
[1981] first suggested that the drifting narrowband features 
were the result of a propagating disturbance which triggers 
narrowband AKR emission at the local cyclotron frequency. 
Several recent papers ascribe AKR fine structure to solitary 
structures such as electron holes [Pottelette et al, 2003], ion 
holes [Mutel et al, 2006], tripolar structures [Pottelette and 
Treumann, 2005] , or to electromagnetic ion-cyclotron waves 
[Menietti et al. , 2006; Xie and Menietti, 2007]. These waves 
propagate along (or obliquely to) field lines, and arc thought 
to modify the electron distribution so as to enhance the CMI 
growth rate. Mutel et al [2006] calculated the effect of ion 
holes on a horseshoe electron and found that the result- 
ing R-X mode CMI emitted power was enhanced by ^^20 
dB. Both the magnitude and the narrowband nature of the 
of the enhancement were consistent with drifting 'striated' 
AKR bursts reported by Menietti et al [1996, 2000]. 

In this paper we investigate three questions stimulated 
by the observational and theoretical investigations of AKR 
summarized above: 

(1) How does the CMI growth rate vary with horse- 
shoe parameters and with propagation direction and am- 
bient density? More particularly, under which conditions 
is the growth a maximum? To answer this, we performed 
a systematic parameter search, varying the electron veloc- 
ity distribution parameters, the ratio of electron plasma to 
cyclotron frequency, and wave-normal angle to the ambient 
magnetic field. 

(2) What is the ratio of growth rates for R-X mode com- 
pared with L-O mode for horseshoe distributions with the 
same range of parameters? Are there parameters for which 
direct L-O mode growth rates reproduce the observed flux 
range of 1% - 10% of R-X mode power? Is there a signifl- 
cant difference in maximum growth rate k-vector orientation 
for LO-mode vs. R-X mode which may be relevant to the 
observations of differing RX/LO ray paths? 

(3) Under what conditions can solitary structures enhance 
(or diminish) the CMI growth rate as a result of the modi- 
fication of the electron distribution function interior to the 
structure? 

Wo address these questions by calculating the growth 
rates and emergent spectral power for a parameterized 
horseshoe electron velocity distribution. The distribution 
function has been designed to closely fit the observed a 
canonical velocity distribution in the Earth's auroral cavity. 
By varying the parameters centered on the best-fit values to 
the observed distribution, as well as the propagation angle, 
we investigate how growth rates vary with each parameter 
and propagation direction, and how the growth rates are 
modified by perturbations . We consider only fundamental 
(flrst harmonic) mode emission and do not consider genera- 
tion of the trapped Z mode, or mode conversion processes. 
We also assume a steady-state process in which the input 
electron beam is constant and hence do consider not the 
evolution of the electron distribution function after interac- 
tion with the radiation field. Finally ,wc assume that the 
fraction of cold electrons in the auroral cavity is very small, 
so that the energetic component dominates the total elec- 
tron density. 

2. CMI growth rates: Parameter search 

2.1. CMI growth rate equations 

Growth rates for the mildly relativistic Doppler-shifted 
cyclotron resonance condition can be written [ Wu and Lee, 
1979; Wu, 1985]: 
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for the right-hand polarized extraordinary (R-X) mode, 
and 
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(3) 



for the left-hand polarized ordinary (L-O) mode. In both 

equations rie is the number density of energetic electrons, 
u)pe is the plasma frequency, / is the normalized velocity 
distribution function (so that the integral of / over velocity 
is the number density), and we have assumed Wr ~ ^ce- We 
also assume nearly perpendicular propagation (fcy <C A:) and 



hence omit the term 



3f 



since it is insignificant compared 



with the df /dv± term. The dimensionless factors fJ-R,L cor- 
rect for the effect of plasma dispersion. This was ignored in 
the original derivation of Wu and Lee [1979], who assumed 
a refractive index n = 1. The values of ij.r,l are discussed 
in Appendix A. 

The delta function is given by the Doppler-shifted gyro- 
resonance condition 



fie 
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(4) 



For V <€. c this equation can be rewritten in the useful 
form 

<5(w||,vx) =<5|^ vl + (vii - Vcf - Vr ^ (5) 



The argument of the delta function is the equation of a 
circle whose center is displaced along the vu axis by 



and whose radius is 



fc||C_ 



(6) 
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where Slj = {ujr — fice) /^ce- The growth rate equations 
(2) and (3) can be recast by converting to polar coordinates, 
resulting in 
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where the integration is performed around the resonant 
semi-circle in {v\\jV±) velocity space. This representation is 
useful in showing that the dependence of the growth rate on 
the velocity distribution is a linear function of the perpen- 



dicular derivative weighted by a simple geometrical factor 
centered at 90° (R-X mode) or 45° and 135° (LO mode). 

2.2. Model velocity distribution function 

We modeled the velocity distribution as a partial annulus 
(horseshoe) with five adjustable parameters: Scale param- 
eter F() (proportional to the number density) , radius {vr), 
characteristic width (a), loss-cone opening angle 6 in the 
upward (anti-Earthward) direction, and horseshoe contrast 
parameter 77. The functional form of the assumed velocity 
distribution function is 



/(^;x,t;||) = Fo|^erfc(^M_ 



77 • exp 



\v\ - Vr 



X g 



(10) 



where the loss cone truncation function g is 

1 if |atan2 (-vx,-v||)| > 9 
otherwise 

where atan2 is the quadrant-dependent inverse tangent 
function. The velocity \v\ is the quadrature sum of v± and 
and erfc is the complementary error function. The veloc- 
ity distribution function is normalized such that the integral 
over velocity space is equal to the electron number density 

f°° 3 
Ue = d vf{v) = 

J —00 

/OO fOO 
dv\\ / dv±v±f [v±,v\i) (11) 
-00 Jo 



The form of the model distribution function, consisting of 
narrow horseshoe population superposed on a flat plateau, 
was chosen to match observed velocity distributions in the 
auroral cavity. Fig. 1(a) shows a pseudo-color plot of a 
model distribution function, while Fig. 1(6) shows a one 
dimensional radial profile comparing an observed distribu- 
tion sampled in the auroral cavity by the FAST spacecraft 
[McFadden et al., 2003] The open diamonds are from Fig. 6 
upper panel, FAST orbit 1804, averaged over 0° - 68° pitch 
angle and converted to velocity units. The closed circles are 
from McFadden et al. Fig. 11 (lowest trace), FAST orbit 
11666. Although the data points are measurements made 
at different epochs, they were both observed while FAST 
was inside the auroral cavity and agree with each other and 
with the model distribution very well. Note that we have 
not fitted the low-energy electron population {E < 100 eV) 
since most are likely photoelectrons and secondaries gener- 
ated by the spacecraft [Delory et al, 1998]. The loss cone 
width angle was fitted to the electron beam depletion evi- 
dent between about 140° — 220° pitch angle [McFadden et 
al, Fig. 6, lower panel]. 

The best-fit model parameter values corresponding to the 
\v\=vr- observed distribution described above are given in the col- 

(8) umn labeled 'FAST match' in Table 1. The corresponding 
electron number density, as given by integrating over the 
distribution function (equation 11), is 0.32 cm^ , which is 
in exact agreement with the measured density reported by 
McFadden et al. 

(9) 2.3. Pcirameter secirch 

'l="r addition to the five parameters associated with the 

model velocity distribution, we also varied the propagation 
direction (parallel wavenumber ratio knjk). We varied all six 
parameters over the ranges listed in Table 1. For each set 
of parameters, the growth rate integrals were evaluated nu- 
merically using Ridder's method for the partial derivatives 
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and Romberg's method for the integrations [Press, et al., 
1995]. In addition, the ahowed frequencies between fice and 
the RX or LO cutoff frequency were divided into 30 intervals 
and the maximum growth rate was recorded. The computa- 
tion of a complete parameter search over the allowed range 
of all variables resulted in computation of growth rates for 
more than 240 million combinations of parameters. In order 
to determine which velocity distributions and ambient con- 
ditions resulted in the largest growth rates, we sorted the 
output datasets separately by RX and LO growth rates. 

2.4. Frequency dependence of growth rate 

Fig. 2(a) shows a pseudo-color plot of the R-X mode 
growth rate integrand (equation 8) for the model distribu- 
tion shown in Fig. 1, along with resonant circles correspond- 
ing to maximum growth rates for perpendicular propagation 
(solid circle), k^^/k = -0.02 (-1.1°), and fcn/fc = +0.02 (1.1°). 
Fig 2(6) shows the corresponding growth rates as a function 
of fractional frequency 



The narrow-band nature of CMI emission is evident, as is 
the strong dependence of growth rate on propagation angle 
(discussed further in section 2.6). The fractional bandwidth, 
which we can define as the full width to half-maximum, is 
about 0.002. This is in close agreement with the model cal- 
culations of [ Yoon and Weatherwax, 1998] who used a similar 
analytic distribution function, and with the numerical sim- 
ulations of [Pritchett et al., 1999]. 

In order to compare with observed fractional bandwidths, 
we must consider the relation between growth rate and emer- 
gent spectral power. The spectral power is very sensitive 
to small changes in growth rate: For a typical c-folding 
amplification factor of 10 (section 2.8.2), a 13% reduction 
from peak growth rate corresponds to a 10 dB intensity de- 
crease. In addition, for non-perpendicular propagation, the 
resonant circle and the velocity distribution form osculat- 
ing circles, so that very small changes in resonant circle ra- 
dius result in large changes in growth rate. Since the radius 
of the resonant circle is a function of fractional bandwidth 
(equation 7), the growth rate integrals are very sensitive to 
Su). Figure 3 shows the spectral power fractional bandwidth 
(defined as full width to -10 dB points) as a function of 
horseshoe energy for a range of propagation angles centered 
on A;|| = 0. The majority of points cluster between 10~^'® 
and 10^^^, but at low energy and non-perpendicular prop- 
agation, some fractional bandwidths are as small as 10"''' . 
These fractional bandwidths encompass most observed nar- 
rowband AKR features, including the isolated striated AKR 
bursts reported by Mutel et al. [2006] who measured frac- 
tional bandwidths as small as 10^^. They are not quite as 
narrow as those reported by Baurnback and Calvert [1987] 
who reported fractional bandwidths as small as 3 x 10~®. 

2.5. RX-mode: Perpendicular propagation 

After inspecting the parameters associated with the high- 
est growth rates, a number of trends became clear. First, 
we examine the dependence on loss-cone opening angle and 
horseshoe width. Figure 4 shows the RX mode growth rate 
as a function of horseshoe width and loss cone opening angle 
for the case fcy = and with other model parameters fixed 
at the 'FAST match' values in Table 1. For loss cone angles 
larger than ~ 90° the growth rate decreases very rapidly, 
indicating that the CMI mechanism requires broad pitch- 
angle scattering of the downward beamed electrons, result- 
ing in a significant population of upward-traveling electrons 
and hence a horseshoe rather than a crescent geometry. 

The horseshoe width parameter, which is a measure of 
the energy spread of the beamed electrons, has a smaller 



effect on the growth rate. At the 'FAST match' coordinates 
(labeled FAST v.d.f.), the width is 0.78 keV. A width change 
of 0.01 (0.6 keV) corresponds to a growth rate change of a 
factor of 2. 

Figure 5 shows R-X mode growth rate as a function of 
horseshoe enhancement factor r? and radius Vr with other pa- 
rameters fixed at 'FAST match' values as above. The right 
ordinate is the electron beam kinetic energy calculated from 

the horseshoe radius. As might bo expected, the growth 
rate increases with increasing beam energy and q. The lat- 
ter parameter characterizes the steepness of the horseshoe 
compared with the plateau component and is proportional 
to the partial derivative df /dv± in the growth rate integral. 

2.6. RX-mode: Non-perpendiculeir propagation 

An illustrative case of non-perpendicular propagation is 
shown in Figure 6. which displays the RX-mode growth 
rate as a function of propagation direction and loss cone 
opening angle, with other parameters fixed at 'FAST match' 
values. The right ordinate displays propagation angle (mea- 
sured from the normal direction) in degrees. The box labeled 
'FAST orbit 1907' shows the estimated loss cone opening an- 
gle (with uncertainty) and k-vcctor propagation direction for 
FAST orbit 1907 observations on 13 Feb 1997 near 18:58:55 
UT. The loss cone opening angle estimate was made using 
the distribution function shown in Fig. 2 of Delory et al. 
[1998] , while the propagation direction and uncertainty was 
taken from Ergun et al. [1998a] for the same time period. 

The growth rate is sharply peaked at = with a 
characteristic FWHM (full-width at half-maximum) ±0.015 
, i.e., a propagation FWHM angle 1.7°. This agrees very 
well with measurements of k-voctor direction observed by 
the FAST electric field sensors. There is also a small but 
significant bias toward the loss cone near 6 — 90°. This is 
a result of the resonant circle intercepting the inner shoul- 
ders of the upward loss cone. While varying all other model 
parameters, we found no circumstances for which the prop- 
agation direction at peak growth deviated significantly from 
perpendicular. 

2.7. LO-mode growth rate 

Figure 7 shows growth rates for RX and LO mode emis- 
sion as a function of horseshoe opening angle using fixed 
'FAST match' for all other model parameters. The Figure 
shows the growth rates in each mode using throe propaga- 
tion directions: fcy/fc = 0°, (solid line), 1° (dashed line), 
and —1° (dotted line). For these parameters, which are typ- 
ical of the auroral cavity, the LO to RX mode growth rate 
ratio is ~ 10~^. We computed the ratio of RX to LO growth 
rates the over the entire range of parameters listed in Table 
1. The ratio was less than 2 • 10 for all RX mode growth 
rates u)i > 10"*. 

Wc can understand this ratio qualitatively by inspection 
of the growth rate expressions (equations 8, 9). At a given 
density, the LO mode growth rate integrand differs only by 
the multiplicative term {vrsin(f>/c)^ = («||/c)^ --^ 0(10^^) 
compared with the RX mode growth integrand. However, 
the remaining integrand heavily weights the line integral to 
regions near v\\ « 0, resulting in a further reduction ~ 10^^. 
Finally, the RX-mode dispersion factor reduces the RX- 
mode growth rate relative to the LO-mode by ~ 10~^ (Ap- 
pendix A). Hence, the expected overall ratio of RX to LO 
growth rates is ~ 10~^. 

2.8. Comparison with observed AKR intensities 
2.8.1. Group speed and convective growth length 

In order to compare calculated growth rates with obser- 
vations of AKR intensities, we must calculate the expected 
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emergent spectral power in each mode. To do this wc need to 
know the physical length over which the CMI maser operates 
and the group speed of the wave. We must also determine 
the background intensity to determine the ratio of emergent 
flux to input flux at a given frequency. 

The temporal growth rate Wi and power gain G (dB) of 
an unsaturated CMI maser are related by 



G{dB) = 10 • logio 



exp 2uj. 



(13) 



where Lc is the convective growth length, i.e., the phys- 
ical length over which the conditions for CMI growth are 
maintained, and Vg is the group speed of the mode. A 
proper calculation the group speed of each mode requires 
solving the relativistic dispersion equation for the assumed 
hot electron velocity distribution including a possible cold 
electron population which may also be present. This calcu- 
lation is beyond the scope of the present paper. However, 
for the R-X mode, we can use the results of Pritchett [1984b] 
who calculated group speeds assuming a DGH distribution 
and a relativistic dispersion relation. He found that the 
group speed varied in the range 0.2c > Vg > 0.1c both for 
a range of DGH parameter 2 < Z < oo and no cold elec- 
trons, and for the case I = 2 and fractional cold electrons 
0.0 < ric/nh < 0.75. Wc adopt an R-X mode group speed 
in the middle of this range, Vg = 0.15c, although the un- 
certainty in this estimate is probably a factor of 2. For the 
LO mode, we assume Vg = c since the radiated frequency is 
much larger than the LO cutoff frequency {u>pe)- 

The convective growth scale is constrained by the phys- 
ical dimensions of the auroral cavity. In the latitudinal di- 
rection the cavity size is ~100 km or less [e.g., Boehrn et al, 
1994]. In the direction along the tangent plane, the auroral 
cavity is much larger. However, the CMI mechanism re- 
quires very similar environmental conditions along the con- 
vective growth path. Probably the most severe limitation 
is the change of magnetic field orientation with respect to 
the propagation vector. For a dipolar field a longitudinal 
shift of 1° at an altitude of 4000 km corresponds to linear 
distance of 180 km. Inspection of Fig. 6 shows that the 
growth rate is reduced by a factor of 2 for angular deviation 
of 1° from maximum growth. Hence we conclude that the 
convective growth scale is unlikely to exceed 200 km in the 
auroral cavity. 
2.8.2. RX Mode Gain 

Omidi and Gumett [1982] estimated that the CMI maser 
required about 10 e-foldings (e^" power gain, or 87 dB) to ac- 
count for the 'observed AKR intensities' , a number that has 
been widely used in the AKR literature subsequently. How- 
ever, more recent in situ obsc^rvation of peak AKR, emission 
electric field intensity in auroral cavities arc now available 
which are significantly higher than previous estimates based 
on measured flux densities from remote spacecraft. Further- 
more, since the sky background intensity and AKR peak in- 
tensities are highly frequency-dependent, we re-examine the 
required maximum gain and its possible frequency depen- 
dence. 

At frequencies between 200 KHz and 800 KHz the sky 
background is dominated by synchrotron emission from the 
galactic plane. The emission peaks near 3 MHz and is 
strongly absorbed by free-free absorption from cold neu- 
tral hydrogen at lower frequencies [Brown, 1973]. The sky 
brightness is steeply rising with an approximate power-law 
form 



Bsky{l^) = -Bo 



(14) 



where Bo = 5 10"^^ W m'^ Hz"^ sr'^ and vo = 360 Khz. 
This dependence is steeper than that expected for free-free 
absorption (exponent 5/2) and may indicate that Razin sup- 
pression by cold plasma in the synchrotron source region is 



significant. Also, since the emission originates in the galac- 
tic plane, the background is not isotropic, but varies by a 
factor ~2 over the sky at a given frequency. The constant 
Bo represents an average over the maximum and minimum 
observed brightness at each frequency [Brown, 1973] . 

In order to compare with observed AKR spectral power, 
we convert the observed sky brightness to squared elec- 
tric field spectral power [Bekefi, 1966; Benson and Fainberg, 
1991] 



Isky(v) — Bakyiv) 

eoC 



(15) 



where Igky is the averaged isotropic spectral power (V^ m~^ 
Hz^^) and eo is the permittivity of free space. For AKR 

radiation measured at FAST orbital altitudes in the auro- 
ral cavity {u ~ 360 KHz), the background flux is 5 10~^* 
(V/mf Hz-i. 

A number of papers have addressed the statistical prop- 
erties of AKR spectral power as measured by remote space- 
craft [e.g.. Green and Gallagher, 1985; Benson and Fainberg, 
1991; Kumamoto et al., 1998; Green et al, 2004]. However, 
in order to estimate the emergent spectral power using flux 
measurements from a remote spacecraft, the angular beam- 
ing pattern of AKR radiation nmst be known. This is cur- 
rently not well established, so we use instead measurements 
of AKR spectral power from in situ spacecraft. The plasma 
wave tracker instrument on the FAST satellite, when op- 
erated in burst waveform mode [Ergun et al., 2001], allows 
direct measurement of AKR intensities within the auroral 
cavity. 

One of the most intense AKR bursts observed by FAST 
had a spectral power 2 10~* (V/m)^ Hz^"*^ [Ergun et al,, 
1998a]. We are not aware of a statistical analysis of AKR 
spectral powers observed by FAST, but Ergun et al. [1998a] 
state that the average spectral power was 'roughly one order 
of magnitude higher than that reported by Viking' [de Fer- 
audy et al., 1987], which implies an average spectral power 
near 10~* (V/m)^ Hz~^. Using these estimates and the 
background spectral power calculated above, the peak CMI 
gain at 360 KHz is 161 dB (18.5 e-foldings) and the aver- 
age gain is 93 dB (10.7 e-foldings). The average is in good 
agreement with earlier estimates, although the peak gain is 
significantly larger. 

Figure 8 is a plot of the calculated RX-mode spectral 
power versus plasma frequency ratio for fcy = 0, group speed 
Vg = 0.15c, and convective growth lengths = 20 km 
(dotted line) and 35 km (dashed line). Wc varied the den- 
sity parameter Fo, but fixed the remaining parameters using 
'FAST match' values (Table 1). The point labeled 'FAST or- 
bit 1907' indicates the observed spectral power and electron 
density (with uncertainty) for the intense burst observed by 
FAST on 13 Feb 1997 at 18:58:55 UT [Ergun et al, 1998a]. 
The point marked 'FAST orbit 11666' is an AKR burst ob- 
served by FAST on 2 Aug 1999 at 23:33:52, and using the 
density measurement reported in McFadden et al. [2003] and 
spectral power of the AKR burst measured from the FAST 
data archive. The observed AKR spectral powers agree very 
well with the model, although the convective growth scales 
should be considered minimum estimates, since any devia- 
tion from normal propagation reduces the growth rate con- 
siderably, necessitating larger convective growth lengths. 
2.8.3. LO mode gain and RX/LO intensity ratios 

We next investigate the question of whether a parameter 
search of calculated LO mode growth rates can yield ratios 
of RX to LO growth rates which are consistent with observed 
values. When LO mode is detected, it is often a few percent 
of the RX mode intensity [Mellott et al, 1984; Benson et 
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al, 1988] although it occasionally reaches 100% [Hanasz et 
al, 2003]. 

Figure 8 shows the LO mode gain using 'FAST match' 
parameters an assumed for group velocity Vg = c and Lc ~ 
200 km (solid line) ), which we have argued (section 2.8,1) is 
the largest plausible convective growth length near 4000 km 
altitude. The gain does not approach observed levels (8-10 
e-foldings) unless the plasma ratio ujpe/i^ce > 0.3. This im- 
plies energetic electron number densities at least 100 times 
larger than observed in the aurora cavity. This suggests that 
LO mode AKR radiation, at least inside auroral cavities, is 
not generated directly by the CMI mechanism. Simulations 
by Pritchett and Strangeway [1985] suggest that a relativis- 
tic effect coupling X-mode (perpendicular) electric fields into 
parallel (0-modo) fields may account for the observed LO to 
RX-mode intensity ratios. The LO mode may also be cre- 
ated from conversion of RX mode emission by polarization 
mismatch in regions of high density gradients [Hayes and 
Melrose, 1986; Louam and Le Queau, 1996b]. 

3. Growth rate modification caused by 
solitary structures and EMIC waves 

We next consider change of CMI growth rate and con- 
sequent gain change caused by modification of an initial 
horseshoe distribution function within solitary structures 
(a.k.a. waves or holes). Solitary structures have recently 
been suggested as the cause of frequency fine structure seen 
in dynamical spectra of AKR radiation [Pottelette et al., 
2001, 2003; Pottelette and Treumann, 2005; Mutel et al., 
2006). We also briefly consider the effect of electromagnetic 
ion-cyclotron (EMIC) waves, which have been associated 
with stimulation of AKR [Menietti et al. , 2006; Xie and 
Menietti, 2007]. 

We assume a simple analytic model for electron and ion 
solitary structures as spherically symmetric potential struc- 
tures with an electric field structure as shown in Figure 9. 
The electron hole accelerates the parallel component of all 
electrons passing the structure, while the ion hole causes a 
reduction in the parallel speed (in both cases, the perpen- 
dicular component is largely unaffected since the gyroradius 
is much smaller than the size of the hole). The parallel ve- 
locity perturbation at the center of the hole can estimated 
by differentiating the energy conservation equation 



SVu 



'meVu 



(16) 



where 5$ is the effective hole potential, and the perturba- 
tion is negative for ion holes and positive for electron holes. 
This local deformation of the distribution function can ei- 
ther enhance or diminish the CMI growth rate depending 
on whether the deformation shifts the perpendicular veloc- 
ity gradient closer to or away from the resonant circle. Even 
if the effect is to enhance the growth rate, the spatial scale of 
the deformation must be sufficiently largo so that the CMI 
gain is significant over the structure. In order to investigate 
the full range of possible conditions, we calculated the RX 
mode growth rates inside and outside the structures, vary- 
ing all model distribution parameters as well as propagation 
direction. 

3.1. Perturbation by Ion holes 

Ion solitary waves (ion holes) are observed in the upward 
current region as symmetric bipolar electric field structures 
with negative potential. They are often seen moving upward 
at several hundred km s^^ in association with upgoing ion 
beams [McFadden et al., 2003]. They have spatial scale sizes 
2-12 km and typical amplitudes £ ~ 10 mV/m, although 



large structures can be several hundred mV per meter. In- 
terior to such large structures, the parallel component of the 
velocity distribution function will be shifted significantly to- 
ward the origin, by up to 1 kcV in energy, corrcsijonding to 
a 5% - 15% decrease in parallel velocity for typical range of 
energetic electron energies. 

We model this perturbation by adding a velocity- 
dependent term to the distribution function 



SVu 



(17) 



where k is a (dimensionless) perturbation parameter. Elec- 
trons with kinetic energies less than the hole potential will 
suffer parallel velocity reversal rather than the decrement 
given in equation 16, as illustrated in Fig. 9b. 

A detailed, self-consistent calculation of the deformed dis- 
tribution function involves solving the Vlasov equation given 
the potential structure of the hole, which is beyond the scope 
of this paper. Instead, we adopt a heuristic approach by ad- 
justing the horseshoe density shoulder factor q to fit the 
observed distribution function within an ion hole. This ap- 
proach is justified by the close agreement between the model 
distribution and the observed distribution within an ion soli- 
tary wave. An excellent example of a distribution function 
perturbed by the passage of a large (700 eV) ion solitary 
wave in shown in Fig. 11 of McFadden et al. [2003]. The 
figure shows a time-sequence of one-dimensional distribution 
function profiles (averaged over 112° — 180° pitch angle) as 
the solitary wave passed by the FAST spacecraft. The se- 
quence clearly shows a 700 eV energy shift toward the origin, 
comparable to the total potential of the wave, as well as a 
3x increase in the peak amplitude. Figure 10 show sampled 
points from this time-series for the unperturbed distribu- 
tion (open squares) and for the maximally shifted distribu- 
tion (filled circles). The dashed line is the best-fit model 
using 'FAST match' parameters (cf. Fig. lb), while the 
solid line was obtained by adjusting k and rj and retaining 
'FAST match' values for all other parameters. We used this 
perturbed model distribution to calculate the growth rate 
inside and outside the solitary wave (see labeled points. Fig. 
5) 

We now assume the effective convective growth length is 
1 /e the overall spatial scale of the solitary wave (Lc ~4 km 
for large structures), and a group speed Vg w 0.15c, as dis- 
cussed in section 2.8. Figure 11(a) shows the resulting RX 

mode power gain as a function of propagation direction for 
the unperturbed state (solid line), and for the perturbed dis- 
tribution within the ion hole (dashed line). The lower plot 
(Fig. 11 (6)) shows the gain difference, which exceeds 40 dB 
over the central part of the gain curve. This is comparable 
with the intensity contrast of the brightest striated AKR 
bursts compared with the background AKR as reported by 
Mutel et al. [2006], which were ascribed to perturbation by 
ion holes. 

3.2. Perturbation by Electron holes 

Electron holes (also called electron solitary structures or 

fast solitary waves) arc isolated regions of positive potential 
with large electric fields, up to 2.5 V/m. They arc nor- 
mally found in the downward current region [Ergun et al., 
1998b]. However, Pottelette and Treumann [2005] argue that 
the presence of 'tripolar' structures in the upward current 
region provides evidence for nested ion and electron holes in 
this region, similar to numerical simulations [Goldman et al, 
, 2003], and that the electron holes are the cause of AKR 
fine structure. 

In order to determine whether electron holes can modify 
AKR growth rates, we calculated the change in growth rate 
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caused by an single electron hole using the same technique as 
the ion hole calculation, viz.,, a spherical distribution (Fig. 
9b) in which incoming electrons are accelerated with a a 
velocity-dependent change given by equation 16, but with 
a positive sign. We determined the growth rate change for 
the entire range of all parameters listed in Table 1. Most 
parameters produced either a negligible difference or a net 
reduction in growth rate, up to 30%. Even in cases of sub- 
stantial reduction in growth rate, the change in emergent 
spectral power was negligible since the overall size of elec- 
tron holes (-^0.5 km, Ergun et al. [1998b]), and hence the 
convective growth scale, is too small to substantially affect 
emergent power. We conclude that electron holes are very 
unlikely to significantly modify AKR emission intensities. 

3.3. Perturbation by EMIC Waves 

Electromagnetic ion-cyclotron (EMIC) waves have many 
properties similar to ion holes: They are frequently detected 
in the auroral cavity, they have peak amplitudes up to ^ 0.5 
keV, and spatial scales (wavelength) 1-10 km [Chaston et al., 
1998]. Similar to ion holes, they also have been observed to 
modulate the ambient electron energy distribution with an 
amplitude comparable to the total EMIC potential [McFad- 
den et al., 1998]. Menietti et al. [2006] and Xie and Mem- 
etti [2007] suggest that EMIC waves could stimulate AKR 
emission by modulation of the electron phase space density. 

From this perspective, an EMIC wave train can be consid- 
ered as a string of ion holes separated by positive potential 
structures of equal amplitude and spatial dimension. Since 
we have shown that ion holes can increase the AKR spectral 
power by up to 40 dB, while positive potentials decrease the 
gain, the net result for each wave crest would be the same 
as an ion hole of similar amplitude, but with a 50% duty 
cycle, i.e., a 3 dB decrease relative to a similar ion hole. 
The aggregate effect will would result in the same gain en- 
hancement but spread over a much larger frequency range. 
Although there are large variations in EMIC wave packet 
sizes and direction, if we assume a wave train has of order 
lO'^ crests, a wavelength of order 1- 10 km, and a k-vcctor 
direction close to the magnetic field direction, the resulting 
z-dircction spatial width is 10^ — 10'^ km. At a mean height 
of 4000 km this results in an effective bandwidth 0.4 - 4 
KHz. This is similar to the bandwidths of many drifting 
AKR bursts seen in AKR dynamical spectra [e.g., Gumett 
et al., 1979; Morozova et al., 2002] 

4. Summary 

We have used a parameterized horseshoe electron veloc- 
ity distribution function to study the dependence of RX and 
LO-mode cyclotron maser instability growth rates on horse- 
shoe geometry, electron density, and propagation direction. 
The distribution function has been closely matched to dis- 
tributions observed by the FAST spacecraft in the auroral 
cavity . We find that growth rates are most strongly af- 
fected by loss-cone opening angle and propagation angle. 
For loss cone opening angles exceeding 90° the growth rates 
drop dramatically, indicating that significant AKR omis- 
sion only occurs in regions for which pitch angle scattering 
of the downward electron beam produces upward traveling 
electrons. The dependence on propagation angle is sharply 
peaked at fcy = with a characteristic full width at half 
maximum of 1.7°. The fractional bandwidth of CMI emis- 
sion, which has a complicated dependence on horseshoe ge- 
ometry and propagation angle, rang ed from lO""*-^ to 10"^ 
but with a concentration near 10~^. The range incorporates 
almost all observed AKR burst bandwidths, including the 
very narrowband striated bursts reported by Mutel et al. 
[2006]. 

We calculated the expected emergent spectral power for 
each mode given the growth rates, frequency-dependent sky 



background, and estimates for the convective growth scales 
and group velocities. The resulting spectral powers for RX- 
mode matched observed FAST observations of AKR power 
with convective growth scales of order 10 km and propaga- 
tion angles 1 — 2°. The LO-mode spectral power is far lower 
than observed values even with a convective growth scale 
as large as 200 km, which we estimate to be the largest 
viable convective growth length in the auroral cavity. By 
performing a systematic parameter search, we find that LO- 
mode spectral power is consistent with observed values only 
for Wpe/^ce > 0.2, i.e., energetic number densities rie > 50 
cm~^ at 4,000 km altitude, far greater than observed values. 
We conclude that LO-mode in the Earth's magnetosphere is 
probably not generated directly, but arises from mode con- 
version of RX-mode radiation. 

Electron velocity distributions are perturbed by passing 
solitary structures (ion and electron holes) and by electro- 
static and electromagnetic waves. We investigated how the 
RX-mode growth rate is affected by an ion hole, by model- 
fitting FAST measurements of the perturbed velocity distri- 
bution caused by an ion hole. We find that large (700 V) 
ion holes can enhance the CMI spectral power by up to 40 
dB, with the maximum gain near fcy = 0. Similarly , we 
suggest that this result can be applied to intense electro- 
magnetic ion-cyclotron (EMIC) waves, which act as a train 
of ion holes. The resulting bandwidth of EMIC-enhanced 
emission is much larger than for an ion hole, typically a few 
kHz. We used a similar analysis for electron holes, search- 
ing all ranges of distribution function parameters. Since all 
growth rates with electron hole perturbations were either 
unchanged or smaller than those without a perturbation, 
we conclude that electron holes cannot be 'elementary radi- 
ation structures' [Pottelette et al., 2001, 2003]. 

5. Appendix A: Correction for relativistic 
dispersion 

The growth rate integrals (equations 2, 3, with jJ,R,L~l) 
were first derived by Wu and Lee [1979], assuming that 
plasma dispersion could bo neglected (refractive index n — 
1). For the ordinary mode and normal propagation, this is a 
good approximation, since the dispersion relation for mildly 
relativistic, zero-temperature plasma only differs from the 
cold-plasma ordinary mode by the relativistic mass [Pritch- 
ett, 1984b] 

u;^ = {kef + ^ (18) 

7 

where 7 is the Lorentz factor of the energetic particles. For 
7 ~ 1 and Wpe << f^ce ~ w the refractive index ul ~ 1, so 
we can set /Uz, = 1. 

The RX-mode dispersion case is more complicated. Lee 
and Wu [1980] generalized the derivation of Wu and Lee to 
include dispersion by a dominant population of cold elec- 
trons. Ormdi, Wu, and Gurnett [1984] showed that the 
RX-mode growth rate equation could be recast, including 
dispersion, by adding a dimensionless multiplicative term in 
front of the resonant integrals 

2r? 

where the function 7?(a;, fc, fice), and the dispersion matrix 
A are given in their Appendix A. They evaluated this term 
assuming that the hot electrons were a small fraction of the 
cold background plasma. They found in a large reduction 
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in the growth rates, up to to two orders of magnitude rela- 
tive to growth rates in the absence of dispersive corrections. 
They concluded that their derived growth rates were too low 
explain observed AKR intensities. 

More recently, spacecraft particle measurements in the 
auroral cavity have shown that the cold plasma assumption 
is invalid: The hot, mildly relativistic component dominates, 
and relatively few (if amy) cold electrons are present [Delory 
et al., 1998]. In addition, in the cold-plasma formalism, RX- 
mode perpendicular waves cannot propagate (since the cut- 
off frequency lies above the gyrofrcquency) , whereas in the 
auroral cavity the most intense radiation is perpendicular. 

In order to solve for the RX-mode growth rate sclf- 
consistently, one must either solve equation (19) using a 
specified distribution function and the relativistic form of 
the dispersion matrix or appeal to numerical simulations 
which incorporate the correct relativistic dispersion rela- 
tions. Pritchett [1986] used such a numerical model, and 
found that the RX-mode growth rates are much larger than 
those computed using a cold-plasma dispersion assumption, 
although still significantly smaller than growth rates com- 
puted in absence of dispersion. 

In this paper we are primarily interested in the depen- 
dence of growth rate on parameterized velocity distribu- 
tions. The conductivity tensor, and hence the dispersion 
matrix, is a function of integrals over the distribution func- 
tion, while the growth rate is most sensitive to the topology 
of the distribution function itself, especially the locations 
of positive velocity gradients with respect to the resonance 
ellipse. Hence to first order we can consider the disper- 
sion correction as a scale factor whose value depends only 
weakly on the model parameters. We have estimated this 
scale factor by comparing the growth rate computed using 
the resonant integral formalism with the numerical simula- 
tions obtained by Pritchett [1984b] using the same distribu- 
tion function, viz., a DGH distribution with parameter range 
2 < Z < 10 and perpendicular propagation. The growth rate 
ratio for the two calculations gave a dispersion correction 
factor hr — 0.14 ± 0.03 where the uncertainty corresponds 
to the ratio variation with parameter I. 

In principle this factor could be different with horseshoe 
distribution functions. However, the excellent agreement of 
RX-mode growth rates computed in this paper using the de- 
rived scale factor fin with numerical simulations using simi- 
lar horseshoe distributions [Pritchett and Strangeway, 1985; 
Pritchett, 1986; Pritchett et al., 1999] indicate that the dis- 
persion factor does not depend strongly on details of the 
velocity distribution geometry. 
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Figure 1. (a) Horseshoe model electron velocity distribution function with parameters providing a best- 
fit to measured distribution function described in McFadden et al. [2003] with parameters listed in Table 
1 (FAST match), (b) One-dimensional cut outside the loss cone, compared with measured distribution 
from McFadden et al, Fig. 6 (open diamonds) and Fig. 11 (lowest trace, filled circles). Note that the 
distribution functions in McFadden et al. have been scaled to velocity units. 
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Figure 2. (a) Pseudo-color representation (blue > 0; red< 0) of R-X mode growth rate integrand 
{sin^(f> X df /dv±) for distribution function shown in Fig. 1. Resonant circles are shown for maximum 
growth rate at propagation directions fc||/fc = —0.02 (—1.1°, dashed line), and fc||/fc — 0.0 (0°, solid line), 
— +0.02 (+1.1°, dotted line), (fe) R-X mode growth rate vs. fractional frequency for propagation 
directions given in (a). 
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Figure 3. Log of spectral power fractional bandwidtii (^cj/fice) vs. electron beam kinetic energy. The 
colors indicate propagation direction, as indicated in right-hand color bar. 
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Figure 4. RX mode growth rate vs. horseshoe width and loss cone opening angle at perpendicular 
propagation, with other parameters fixed at 'FAST match' values list in Table 1. Contours lines increase 
by a factor of two starting from 0.0015. The point labeled 'FAST v.d.f.' corresponds to the best-fit 
parameters matching the velocity distribution function described in McFadden et al. [2003], FAST orbit 
1804. 
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Figure 5. Same as Fig. 4 but for horseshoe radius Vr/c and enhancement factor 77. Contour lines 
increase by 10% starting at 0.0488. The point labeled 'FAST v.d.f. (unperturbed)' corresponds to the 
'FAST match' parameters in Table 1. The point 'Perturbed (ion hole)' corresponds to perturbation of 
the velocity distribution by a 700 eV ion solitary wave (Section 3 
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Figure 6. Same as Fig. 4 but for loss cone opening angle 9 and propagation direction fc||/fc. The dotted 
box labeled 'FAST orbit 1907' indicates the observed AKR propagation direction with uncertainties 
[Ergun et al, 1998a] and loss cone opening angle estimated from 2-dimensional distribution function 
[Delory et al, 1998], FAST orbit 1907, 13 Feb 1997 near 18:58:55 UT. 
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Figure 7. Growth rate for R-X mode (filled circles) and L-O mode (open circles) versus loss cone 
opening angle with other model parameters fixed at 'FAST match' values. For each mode, the solid 
lines are for perpendicular propagation (0 = 0°) while the dashed and dotted lines are for ^ = —1° and 
<f) = +1° respectively. 
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Figure 8. Computed RX mode spectral power (left axis) and equivalent c-folding field amplification 
(right axis) of CMI maser as a function of plasma frequency ratio (bottom axis) and number density at 
4000 km {fee = 360 KHz) using model parameters listed in Table 1 (FAST match column), fey = 0, group 
velocity Vg = 0.15c and convective growth length Lc = 35 km (dashed line) and Lc = 20 km (dotted 
hne). The points labeled FAST orbit 11666 and FAST 1907 are measured AKR spectral powers and 
electron densities (see text for details). The LO mode spectral power, using Vg = c and Lc =200 km 
(solid line), does not approach observed intensities until the electron number density is more than 100 x 
the observed mean auroral cavity density. 
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Figure 9. Effect of ion and electron holes on velocity distribution function. ; (a) Model E-field vectors 
for spherical ion hole (top) and electron hole (bottom); (6) Velocity profile versus distance from hole 
center for incident electrons with initial with speeds of 30,000 km s~^ (red), 20,000 km s~^ (blue), and 
10,000 km s~^ (green); (c) perturbed velocity distributions modified by ion hole and electron holes with 
a total potential of 0.7 keV. 
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Figure 10. One-dimensional radial profile of model distribution function using 'FAST match' parameters 
(dashed line). Open squares are FAST observations scaled from one-dimensional distribution function 
plotted in McFadden et al. [2003] Fig. 11 (lowest trace). The perturbed distribution profile (solid line) 
corresponds to velocity best-fit perturbation parameter k = 0.014 and rj = 4.4, The solid circles are data 
points scaled from the middle trace of McFadden et al. when the perturbation from a passing ion solitary 
structure is maximal. 
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Figure 11. (a) R-X mode gain (left axis) and corresponding c-foldings (right axis) versus propagation 
direction for horseshoe velocity distribution function with 'FAST match' parameters list in Table 1 
assuming a 40 km convective growth path and group velocity Vg = 0.15c. The solid line is for the 
unperturbed distribution, while the dashed line is for a 0.7 keV ion hole and the dotted line is a 0.7 
keV electron hole. (6) As above, but gain difference between unperturbed distribution and perturbation 
caused by ion hole. 
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Table 1. Model Parameters and Search Ranges 



Model Parameter 


Symbol 


FAST match^ 


FAST value*" 


Range 


Range Notes 


Density normalization 


Fo 


lQ-18 


0.32 cm"-'* 


1Q-19_ 10-17 


Units: s'^m"*^ 


Horseshoe radius 


Vr 


0.118c 


3.5 keV 


O.lOc- 0.20c 


2.5 < KE < 10.2 keV 


Horseshoe width 


a 


0.013c 


0.78 keV 


0.01c -0.02c 




Loss-cone opening angle 





40° 




0° - 180° 


Full width is 2 X 61 


Horseshoe shoulder 


V 


1.6 




1.0 - 3.0 




Parallel wavenumber ratio 






0"= 


-0.1- +0.1 




Perturbation parameter 


Energy perturbation factor 


K 


0.014 






Eqn. 17 


Shoulder perturbation factor 


v 


4.4 






see FAST match Fig. 10 



* Parameters which best match observed FAST distribution function shown in McFadden et al. [2003] 
^ Physical value corresponding to FAST model parameter. 

" Measured AKR k- vector direction , 0.06 but consistent with fcy = 0, Ergun et al. [1998a], FAST orbit 1907. 
Measured ion solitary wave amplitude McFadden et al. [2003], FAST orbit 11666. 



